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ABSTRACT 

We present high angular resolution (0''7) observations made with the Very 
Large Array (VLA) of the radio recombination line (RRL) H53a and radio con- 
tinuum emission at 43 GHz from the ultraluminous infrared galaxy (ULIRG) 
Arp 220. The 43 GHz continuum emission shows a compact structure (~ 2") 
with two peaks separated by ~ 1", the East (E) and West (W) components, that 
correspond to each galactic nucleus of the merger. The spectral indices for both 
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the E and W components, using radio continuum images at 8.3 and 43 GHz are 
typical of synchrotron emission (a ~ —1.0). Our 43 GHz continuum and H53a 
line observations confirm the flux densities predicted by the models proposed by 



Anantharamaiah et al.l (120001 ). This agreement with the models implies the pres- 
ence of high-density (~ 10^ cm~^) compact HIT regions (~ 0.1 pc) in Arp 220. 
The integrated H53a line emission is stronger toward the non-thermal radio con- 
tinuum peaks, which are also coincident with the peaks of molecular emission of 
the H2CO. The coincidence between the integrated H53a and the H2CO maser 
line emission suggests that the recent star forming regions, traced by the high den- 
sity gas, are located mainly in regions that are close to the two radio continuum 
peaks. A velocity gradient of ~ 0.30 km s~^ pc~^ in the H53q; RRL is observed 
toward the E component and a second velocity gradient of ~ 0.15 km s~^ pc~^ is 
detected toward the W component. The orientations of these velocity gradients 
are in agreement with previous CO, HI and OH observations. The kinematics 
of the high-density ionized gas traced by the H53q; line are consistent with two 
counter rotating disks as suggested by the CO and HI observations. 

Subject headings: galaxies: general — galaxies: individual (Arp 220) — galaxies: 
starburst — radio continuum: galaxies — radio lines: galaxies 



INTRODUCTION 



Arp 220 (IC 4553/4; UGC 9913; IRAS 15327+2340) has an infrared luminosity of 
~ 1.5 X 10^^ L0, as determined from IRAS 25 — 100 /im observations (Soifer et al. 1984), and 
therefore is classified as an ULIRG. At a distance of ~ 70 Mpc (H^ = 75 km s^^ Mpc~^), 
Arp 220 is characterized in the optical by tidal tails interpreted as the signature of a merging 
process (Surace, Sanders fc Evans 2000). The large visu al extinction toward Arp 220, which 
could be Av > 100 mag (iGonzalez- Alfonso et al.ll2004l ). precludes the observations of the 
nuclear regions at optical and even at infrared frequencies. Radio frequency observations 
do not suffer from dust extinction, providing a powerful tool to investigate the kinematics, 
morphology and physical properties of the gas in the dust-obscured regions of Arp 220. Radio 
observa tions reveal a double-nucleus system separ ated by ~1" (~ 350 pc) as observed in CO 
and HI (IMundell et al.ll200ll : ISakamoto et al.lll999l ). The existence of multiple luminous radio 
supernov ae (RSN) in the nu clear regions of Arp 220 suggests that it is mainly starburst- 
powered (ISmith et al.lll998l ). However, compact OH maser and X-rays emission suggest 



the presence of AGN activity in the W nucleus (INorris et al.l Il985l : iLonsdale et al.l Il998 
Clements et alll2002[ l. 
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High angular resolution (0''5) CO observations carried out by ISakamoto et al.l (119991 ) 
suggest that the two nuclei are counter-rotating with respect to each other and are em- 
bedded in a kiloparsec-size gas disk, which rotates around their dynamical center. The 
counter-r otation is also observ ed in higher angular resolution {0'.'22) observations of HI ab- 
sorption jMundell et al.lboOlh . which show that the rotation is not coplanar and proposed 
a model in which the two nuclei are in the final state of merging. Counter-rotation may 
provide the mechanism to get rid of angular momentum, a condition for the two nuclei to 
merge. The senses of the ve locity gradients, as determined from OH maser emission obser- 
vations ( Rovilos et al. 12003), are consistent with CO mo lecular emission and HI absorption 



observations (ISakamoto et al.lll999l : iMundell et al.ll200ll ) 



Arp 220 is the most distan t syst e m from which radio recomb ination lines (RRLs) 
have been observed ( Zhao et al. 1996 : Anantharamaiah et al. 200ol). VLA o bserva tions 

of the centimeter wav e lengt h RRL H92q; have been reported by IZhao et al.l (|l996r) and 

Anantharamaiah et al.l (120001 ) with angular resolutions of 4" and 1", respectively. lAnantharamaiah et al 

using a single dish telescope (IRAM 30m), detected the emission of the millimeter 
wavelength RRLs H42a;, H 40q;, and H31q: and with the V LA established upper limits for the 
RRLs H167a and H165a. lAnantharamaiah et al.l (120001 ) constructed a model using three 
density components of ionized gas in order to reproduce the observed RRL and radio contin- 
uum intensities, obtaining a star formation rate (SFR) in Arp 220 of ~ 240 yr~^. The 
existence of a high-density component of ionized gas, as deduced from these models, has 
been interpreted as evidence for recent star formation in Arp 220. In particular, the RRL 
H53a can be used as a tracer of thermal high-density gas (~ 10^ cm~^). The star formation 
in this system may be supported by multiple starbursts with a high star formation rate. 



In this paper we present VLA observations of Arp 220 in the millimeter regime of the 
RRL H53a (7 mm) as well as the radio continuum at 43 GHz with subarcsecond (0''7) 
angular resolution. The observations are described in §2. In Section 3 the results obtained 
from the 43 GHz radio continuum and the RRL H53a are presented. In § 4 the implications 
of these results are discussed. The conclusions are presented in § 5. 



2. VLA Observations of the H53a line. 

The RRL H53a (43 GHz) was observed with the VLA in the C configuration on October 
13 and 25, 2002. The fiux density scales were determined using observations of J1331-I-305 
(1.49 Jy). The fast-switching mode was used in order to track the phase variations induced 
by the troposphere using cycles of 230 s on Arp 220 and 40 s on the phase reference source 
J1540+147 (~ 1 Jy). The bandpass response of the instrument is frequency dependent 
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and observations of a bandpass calibrator were required for each frequency window. This 
bandpass response was corrected using observations of J1229+020 (17.5 ± 0.5 Jy). Five 
frequency windows, each of these with 15 spectral channels (~ 22 km s~^ each), were used 
to observe the broad RRL. The total velocity coverage is ~1000 km s~^. The different 
frequency windows were centered at 42114.9, 42135.1, 42164.9, 42185.1 and 42214.9 MHz. 
The on-source integration time was ~2 hrs for each frequency window. To avoid possible 
phase decorrelation at 7 mm, the calibration of the data was performed correcting for the 
phases in a first step and subsequently correcting for both amplitude and phase. The spectral 
line data were further calibrated by applying the solutions obtained from the self-calibration 
performed on the continuum channel (which contains the average of the central 75% of the 
bandpass) of each frequency window. The spectral line observations of the second day were 
not combined with those of the first day since the variable content of water vapor in the 
troposphere caused decorrelation of phases in the second and third frequency windows. A 
critical step in the reduction of the data is to determine the bandpass shape by normalizing 
the bandpass using the channel data. Five line images, one for each frequency window 
of the first day, were produced using a weighting scheme intermediate between natural and 
uniform to obtain a circular beam of 0''7. The five line cubes were regridded in frequency and 
combined into a single line cube in the GIPSY environment . The method used to combine 
different windows is described by lRodriguez-Rico et al.l (120041 ). The data were then Hanning- 
smoothed to minimize the Gibbs effect and a final velocity resolution of ~ 44 km s~^ was 
achieved. The parameters of the observations are given in Table 1. The resulting line cube 
was further processed using IMLIN from AIPS to subtract the continuum emission. A first 
order polynomial fit was used for the line-free channels, two on the low-velocity edge and 
five on the high velocity-edge of the bandpass. Previous Ar p 220 observations of the RRL 
H42a (lAnaritharamaiah et al.l 120001) . CO (2-1) emission line (jSakamoto et al.l 119991 ) and HI 
absorption (IMundell et al.l l200ll ) , show that the full width at zero intensity expected for 
the H53a line is in the range from ~ 5000 to 5700 km s^^. The velocity range covered 
by the H53a line observations is 5100 — 6100 km s~^, consequently loosing the continuum 
that corresponds to the blueshifted side of the line. We performed a test using the line- 
free channels on the blueshifted side of the line as a baseline. We have determined that if 
weak H53a line emission with intensity < 3 mJy (3cr) is present in the velocity range of 
5000 — 5100 km s~^ the full width at half maximum (FWHM) is not affected. The fit for 
the H53q; line is in agreement with the fit obtained using the continuum for the low velocity 
edge of the bandpass. The H53 q; (43 GHz) data have been c ompared with the RRL H92a 
(8.3 GHz) observations taken by lAnantharamaiah et al.l (120001 ) with the VLA in the B array. 
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3. RESULTS 

3.1. Radio continuum. 

Figure [T] shows the continuum image of Arp 220 at 43 GHz superposed on the continuum 
image at 8.3 GHz, both images with an angular resolution of 0'.'7. In this image, the two 
continuum peaks (E and W) are clearly resolved at both frequencies. Using the task IMFIT 
in AIPS we estimated the properties of these two components. The physical parameters at 
8.3 and 43 GHz were obtained integrating over equivalent areas in order to estimate the 
spectral index a {S^, oc z/") of each component. 

Table 2 lists the physical parameters of the continuum components E and W determined 
at 43 GHz. The total measured 43 GHz continuum flux density of Arp 220 is 44 ± 4 mJy. 
The E and W components have continuum flux densities of ~ 17 and ~ 21 mJy, respectively. 
There is also a contribution of ~ 6 mJy to the total flux density from extended emission 
that was determined using the task IRING in the AIPS environment by integrating the flux 
density over an angular scale of 3" to 8". This total continuum flux density at 43 GHz is in 
agreement with the value obtained by interpolating from previous observations at 32 GHz and 
97 GHz (Downes & Solomon 1998, Anantharamaiah et al. 2000). The continuum physical 
parameters at 8.3 GHz of both the E and W components in Arp 220 are listed in Table 3. 
The positions of each continuum peak are coincident at 8.3 and 43 GHz. The two radio 
nuclei of Arp 220 are separated by 0''9 ± 0''1 at a position angle of 95° (E with respect to 
N). The spectral index between these two frequencies for the E component is —0.96 ± 0.08 
and for the W component is —0.97 ±0.07. These spectral indices are typical of non-thermal 
emission. 



3.2. Radio Recombination Line H53q; 



We detected H53a RRL emission (> 3 mJy) toward both components, E and W, of 
Arp 220 in the velocity range 5370 to 5845 km s^^. The velocity range over which the H53a 
l ine emission is observed agr ees with observations of the RRLs H40q:, H 42a, and H92a 
(lAnantharamaiah et al.ll2000l ) and the C0(2-l) line (ISakamoto et al.lll999l ). Figure [2] shows 
the VLA spectral channel images of H53a line emission observed over the central ~ 500 pc 
of Arp 220. Figure [3] shows the velocity-integrated line emission (moment 0) of the H53a line 
superposed on the moment of the H92a line. The peaks in the moment of the H53a line 
emission are nearly coincident (< 0''5) with the two radio continuum nuclei. The moment 
of the H92q! line emission (shown in gray scale) exhibits two peaks that are close to the 
continuum peaks and are coincident with the H53a! peaks within 0''2. 
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The total integrated spectrum over the region with detectable H53a line emission is 
shown in Figure H] along with the E and W spectra obtained by integrating over each of 
these components. A single Gaussian fit was carried out for each line profile. The results 
from the H53a line fits are listed in Table 4. The peak line flux density in column (2), The 
FWHM of the line is given in column (3), the heliocentric velocity of the line center in column 
(4) and the integrated line emission in column (5). Based on the fit parameters for the E and 
W components, the peak H53q! line flux density for the W component is approximately 30% 
stronger than the E component. The H53a line widths and central velocities are similar for 
the E and W components. The H53a and H92a line widths and central velocities determined 
for the total integrated line emission over Arp 220 are in good agreement with each other. 
Figure [5] shows the re s ults f rom radio continuum and RRL emission models taken from 



Anantharamaiah et al.l (120001 ). The measurements for the total integrated 43 GHz continuum 



and H53a line emission are shown in Figure andEb, respectively. 

Figure [6] shows the H53a line velocity field (moment 1) for the ionized gas in Arp 220. 
Based on this H53a velocity field, two velocity gradients are observed toward the E and W 
components. The velocity gradient on the E component is ~ 0.30±0.10 km s~^ pc~^, oriented 
along a P.A. of ~ 30°. The velocity gradient on the W component has been marginally 
detected (~ 0.15 ± 0.10 km s~^ pc~^) along a P.A. of ~ 260°. Due to the lower velocity 
resolution of the H92q; line (230 km s~^), compared to the velocity resolution of the H53a 
line (44 km s~^), a detailed comparison of the H53q; and the H92q; velocity fields is not 
possible. 



4. DISCUSSION 



The measured separation between the two continuum peaks (0''9±0''1 ~ 320 pc) and the 
P.A. of their separation (95°) at 4 3 GHz are in good agreemen t with previous measure ments 
(~ 350 pc) at 5, 15 and 23 GHz jNorrislll988h and 4.8 GHz jBaan fc Haschicklll995h . The 
peak positions at 8.3 and 43 GHz agree at the level of 0''1. The spectral index for Arp 220, 
determined using the radio continuum flux density of ~ 170 mJy measured at 8.3 GHz and 
the measurement at 43 GHz of ~ 44 mJy, is a ~ —1.0. On the other hand, the spectral 
index determined using the radio continuum flux density of ~ 49 mJy measured at 32 GHz 
(Baan, Gusten & Haschick 1986) and the measu rement at 43 GHz (this work) of ~ 44 mJy, is 
a ~ 0.36 (±0.05). lAnantharamaiah et al.l (120001 ) estimated that the thermal radio continuum 
flux density at 5 GHz is ~ 32 mJy. Using a spectral index value a = —0.1 for the thermal 
free-free emission the expected thermal flux density at 43 GHz is ~ 25 mJy, which is ~ 60% 
of the total measured continuum flux density at 43 GHz (see Figure 5). The flattening of the 
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spectral index, a value of a closer to zero when the 32 and 43 GHz measurements are used, 
indicates that the ratio between the thermal emission and non-thermal emission is larger at 
43 GHz than at 8.3 GHz. 



Anantharamaiah et al.l (120001 ) modeled the er nission of radio continuum an d RRL ob- 
servations at different frequencies. In these models lAnantharamaiah et al.l (120001 ) used three 
thermally ionized gas components (Al, Bl and D) and non-thermal emission. Compo- 
nents Al and D are characterized by electron densities of 10^ cm~^, while component Bl 
has a larger electron density of 2.5 x 10^ cm~^. All p ararn eters for each of these com- 
ponents are listed in Table 9 of lAnantharamaiah et al.l (120001 ). Figure [5] shows the model 
results from lAnantharamaiah et al.l (120001 ) including the measurements of the H53a; line 
(~ 510 X 10~^^ W m~^ or 3.7 Jy km s~^) and the continuum flux density at 43 GHz 
(~ 44 mJy). The H53a line and 43 GHz continuum fl ux densities are in agreement with the 
values predicted in the lAnantharamaiah et al.l (120001 ) models, conflrming the presence of a 



high-density ionized gas component (~ 10^ cm~^) composed of compact HII regions with 
~ 0.1 pc diameter. According to these models the contribution from the dense gas to the 
H92a line emission is negligible and the contribution from the lower density gas to the H53a 
line is < 10%. 



On the other hand, the models of lAnantharamaiah et al.l (120001 ) indicate that ~ 50% 
of the H53q; line emission in the high-density ionized gas is due to internal stimulated line 
emission. The spatial distribution of the H53q! and the H92q; lines is similar toward the 
E component as can be seen in Figure [31 However, a remarkable difference is observed 
toward the W component where the H92a line emission is more extended than the H53a 
line emission. Near the continuum peaks, the spatial distribution of the H53a and H92a 
RRLs emission suggests that the recent star formation is taking place mostly in two regions. 
These two regions of line emission are slightly displaced from the respective continuum peaks: 
the W line peak is 0!'2 ± OCl to the S of its associated continuum peak and the E line peak 
is 0!'3 ± 0''1 to the N of its associated continuum peak. According to the integrated H53a 
line emission (see Figure [3]), the bulk of high-density ionized gas is concentrated in these 
two regions. This result is supported by the spatial correlation between the two velocity- 
integrated H53a line emission pe aks and the peaks of the formaldehyde (H2CO) molecular 
emission (iBaan fc Haschicklll995l ). 



The models of lAnantharamaiah et al.l (120001 ) derived a SFR in Arp 220 of 240 yr~^, 
assuming a mass range of 1 — 100 in the Miller-Scalo initial mass function (IMF). 
According to these models, the SFR could be as low as 90 yr~^ if the upper mass 
limit in the IMF is reduced to 60 M© and the Salpeter IMF is used. Given that the radio 
continuum and the velocity-integrated line flux density are in agreement with the values 
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expected from the models of lAnantharamaiah et al.l (120001 ). the value for the SFR could be 



in the range of 90 to 240 M0 yr-i and the mass of high-density ionized gas (~ 10^ cm ^) is 
between 10^ and 10"^ Mq. 

The orientations of the H53q; velocity gradients in Arp 220, with P.A.~ 30° in the E 
component and P.A.^ 260° in the W corn ponent agree with the CO and HI observations 
(ISakamoto et al.lll999l : iMundell et al.ll200ll ). Based on t he P.A. of these t wo ve locity gradi- 
ents, the CO (2-1) line observations were interpreted by ISakamoto et al.l (119991 ) as evidence 
of two counter-rotating disks in Arp 220. An alternative model that consis ts of a warped gas 
disk th at resulted from a merger of two spiral galaxies was also proposed by lEckart &: Downes 
(I2OOII ). The angular resolution achieved in the H53a line observations (0''7) is insufficient 
to discern between these two models. The velocity structure of the ionized gas based on 
the H53a data show two velocity gradients in the two components of Arp 220. On the E 
component the H53a velocity gradient is ~ 0.3 km s~^ The HI velocity gradient 

on the E component is 1.01 ± 0.02 km s~^ pc~^ jMundell et al. 2001) and the OH velocity 

' 2003h. 



gradient is 0.32 ± 0.03 km s pc (iRovilos et al 



The velocity gradient on the E component implies a virial mass of ~ 8 x 10^ (sin^'^i) Mq 
in a disk of radius ~180 pc (z is the inclination). T he total mass calcul ated from OH obser- 
vations, for a disk of radius ~ 80 pc, is ~ 10'' Mq (IRovilos et al.ll2003l ). The value derived 
from the OH observations is consitent with the H53q; line observations when the different 



sizes are taken into account. Based on HI and CO (ISakamoto et al.l Il999l : iMundell et al. 



2OOII ) observations, the estimates of the total mass are ~ 10^ Mq. The total mass estimated 



using CO and HI observations is two orders of magnitude larger than the total mass calcu- 
lated from these H53a line observations. The smaller estimates obtained from the H53a line 
and OH observations are explained if the molecular CO and neutral HI gas are distributed 
over a larger region compared to the ionized gas. In the direction of the W component, the 
H53a vel ocity gradient is a fa c tor of 10 less than o bserved in HI absorption and OH maser 
emission (IMundell et al.ll200ll : IRovilos et al.l l2003l ) . However, the orientation and sense of 
the velocity g radient in H53q; are consistent with those determined from the CO, HI and OH 
observations (ISakamoto et al.l Il999l : IMundell et al.l 1200 ll : IRovilos et al.ll2003l ). 



CONCLUSIONS. 



We have observed the RRL H53a and radio continuum at 43 GHz toward the ULIRG 
Arp 220 with high angular resolution (0''7) using the VLA. The 43 GHz radio continuum 
and the RRL H 53a have been compared wit h observations of the 8.3 GHz radio continuum 
and H92a line (lAnantharamaiah et al.ll2000l ). 



- 9 - 



• The total 43 GHz radio continuum flux density of Arp 220 is 44 ± 4 m Jy. The morpho- 
logical characteristics observed in the radio continuum at 43 GHz agree with previous 
radio observations made in the range from 1.4 to 23 GHz. In the radio continuum at 
43 GHz, Arp 220 exhibits a double-nucleus system that has been clearly resolved with 
a separation of the two radio continuum peaks of ~ 1" at P. A. of 95°. 

• We have determined the spectral indices for both the E and W components, using radio 
continuum flux densities at 8.3 and 43 GHz. Both components have spectral indices 
typical of synchrotron emission {a ~ —0.9), as expected from extrapolation from the 
centimeter wavelength range. 

Based on 43 GHz VLA dat a, we confirm the predictions of the models made by 
Anantharamaiah et al.l (120001 ) for RRLs. The integrated H53q! line flux density is 
about a factor of ~ 50 times la rger than the integrated H92 q: line flux density, in 
agreement with the predictions of I Anantharamaiah et al.l (120001 ). Thus, the H53a line 



traces the high-density (~ 10^ cm ^) compact HII regions (~ 0.1 pc) in Arp 220. 

The spatial distributions of the H53a and the H92a line emission are similar in the 
direction of the E component. On the W component the low density (~ 10'^ cm~^) 
ionized gas component is more extended than the high density (~ 10^ cm~^) compo- 
nent. 



The kinematic and spatial distribution behavior as observed i n the RRL H53q; is in 



agreement with resu lts reported from CO and HI observations (ISakamoto et al.lll999 



Mundell et al.ll200ll ). supporting the counter-rotation of two disks in Arp 220 at smaller 
scales. 



The original version was published in the Astrophysical Journal, volume 633, page 198 in 
year 2005. This new version is submitted with the corrected velocity scale. In Figures 2, 4 and 
6 the values for radial velocities were incorrect. All velocities were increased by 245 km s~^ to 
correct values to heliocentric values (optical definition). To obtain the velocities using radio 
definition, add 145 km s^^. Thus the central velocities in Table 4 should be 5495 km s~^ 
(optical definition) or 5395 km s~^ (radio definition). This velocity shift was introduced 
in the frequency reprojection process. The velocity definition in Table 4 is the optical 
heliocentric definition, as is used by K. R. Anantharamaiah et al. (ApJ, 537, 613 [2000]). 
All other quantities in the paper remain unchanged, e.g. the radio recombination line and 
radio continuum integrated properties, as well as distributions. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. CR and YG 
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Table 1. OBSERVATIONAL PARAMETERS FOR Arp 220 



Parameter H53q; Line 



Right ascension (J2000) 15 34 57.28 

Declination (J2000) 23 30 11.9 

Angular resolution 0''7 

Total observing duration (hr) 13 

Bandwidth (MHz) 150 

Number of spectral channels 45 

Center Y Hei (km s-^) 5500 

Velocity coverage (km s"-*^) 1000 

Velocity resolution (km s~^) 44 

Amplitude calibrator J1331+305 

Phase calibrator J1540+147 

Bandpass calibrator J1229+020 

RMS line noise per channel (mJy/beam) 1 
RMS, continuum (mJy/beam) 0.3 
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Table 2. RESULTS FROM GAUSSIAN FITTING TO THE CONTINUUM EMISSION AT 43 GHZ 

OF ARP 220. 



Feature 


RA (J2000) 


DEC (J2000) 


Size^ , RA. 


Sc (mJy) 


Arp 220 E 


15''34™57f28 ± OfOl 


23°30'11"3±0"1 


0"5 X 0"3, 89° 


17± 2 


Arp 220 W 


15''34'"57f22 ± OfOl 


23°30'llM±0n 


O'fi X 0'.'2, 114° 


21 ±2 



Deconvolved angular size. 
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Table 3. RESULTS FROM GAUSSIAN FITTING TO THE CONTINUUM EMISSION AT 8.3 GHZ 

OF ARP 220^ . 



Feature 


RA (J2000) 


DEC (J2000) Size*^ , RA. 


Sc (mJy)'= 


Arp 220 E 


15''34™57f29 ± OfOl 


23°30'11"3 ± 0"1 0"5 X 0('4, 88° 


77 ± 1 


Arp 220 W 


15''34'"57f22 ± OfOl 


23°30'11^'5 ± 0"1 0"5 X 0"3, 107° 


89 ± 1 



'^Observations at 8.3 GHz made by lAnantharamaiah et al.l (120001 ) . 
'^Deconvolved angular size. 

'^Continuum flux densities were measured in the same area as the 43 GHz continuum 
flux densities listed in Table 2. 
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Table 4. RESULTS FROM GAUSSIAN FITTING TO THE H53a LINE EMISSION OF ARP 220. 



Feature 




A YpWHM 




1.07 (Sp ^Vfwhm ) 




(mJy) 


(km s~^) 


(km s~^) 


(W m-2 X 10-23 ) 


Arp 220 E .... 


6±1 


235 ± 20 


5485 ± 10 


220 ± 40 


Arp 220 W ... 


8±1 


265 ± 30 


5515 ± 10 


320 ± 50 


Arp 220 E+W 


16 ±2 


230 ± 20 


5495 ± 10 


560 ± 80 
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Fig. 1. — Radio continuum image of Arp 220 at 43 GHz (contours) made using VLA obser- 
vations in the C array superposed on the continuum image at 8.3 GHz (gray scale) using 
the VLA in the B array, both images made with angular resolution of 0'.'7. The contour 
levels (43 GHz) are drawn at -3, 3, 5, 7, 9, 18, and 36 times the rms noise 0.3 mJy beam~^. 
The gray scale (8.3 GHz) ranges from 0.15 to 70 mJy beam"-*^. The crosses mark the peak 
positions (given in Table 2) of the 43 GHz continuum. The size of the crosses are 0''5, which 
is about five times the error position of the radio continuum peaks. 
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Fig. 2. — H53a continuum subtracted channel images from Arp 220 labeled according to 
the central heliocentric velocity of each channel. Contours are -3, 3, 4, 5, 6, 7, 8, 9, 10 and 
11 times 1 mJy beam^^, the rms noise. The crosses show the position of the two 43 GHz 
continuum peaks. The synthesized beam is 0''7 FWHM. 
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Fig. 3. — Integrated H53« line emission (moment 0) from Arp 220 in contours with integrated 
H92q! line emission (moment 0) superposed in gray scale. Contour levels are 10, 30, 50, 
70, and 90% of the peak (1.3 Jy beam~^ km s~^). The gray scale covers the range — 
200 Jy beam~^ km s~^. The crosses show the position of the two 43 GHz continuum peaks. 
The synthesized beam for both RRLs is 0(7 FWHM. 



-19- 



^2 



1:1 
a 



Eh 



'tn 

d 

V 

a 

E 



J. ^ 

'tn 

d 

V 

a 

E 



1 

A ^ 1 
/ \^ 1 
/ A / 

/ / V 

h 


J, 1 1 

v 


Arp 220 
H53a 


T 

fi 


\\ 




\ 1 
\ j 




'T~\~/ 


\' 






5445 


5645 5845 

Velocity (km/s) 


6045 6245 


1 


1 1 


1 1 

E 


/ / 
- / 

1/ 


\\ 


_ 


'/ 

/ ^- 






1 i ^ 
/ \ 
-I \y 








\/ \ / 




5445 


5645 5845 

Velocity (km/s) 


6045 6245 


1 J 

/\ ^ 1 

h 


^\ 1 1 
^\ 

\\ 
\ \ 
\\ 

\ 


1 1 

w 


,/ \ /■ 
\ 1 


\ ' i.y / f- - /V- . 

v.^ ,' ,' \ / \ 

\! \ J \ 




5445 


5645 5845 


6045 6245 



Velocity (km/s) 



Fig. 4. — H53q! recombination line profiles in Arp 220 observed using the VLA integrated 
over the region shown in contours in Figure 1 (top), over the E component (middle) and 
over the W component (bottom). The dashed line shows the resulting Gaussian fit to the 
data and the dashed-dotted hne indicates the residuals from the Gaussian fit. 
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Fig. 5. — Continuum and RRL model results for three ionized gas components (Al, A2, and 
D) in Arp 220 by Anantharamaiah et al. (2000). The recombination line and continuum 
flux densities are shown in (a) and (b), respectively. The contribution from the different 

components are shown and the parameters for each component are given in Table 9 of 
Anantharamaiah et al. (2000). The observed line data points are listed in Table 2 and 4 of 
Anantharamaiah et al. (2000). The filled triangles indicate the velocity integrated H53a line 
fiux density in (a) and the 43 GHz continuum fiux density in (b) obtained in the current 
observations. The size of the triangles represent the la error for the measured values. 
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Fig. 6. — Velocity field of Arp 220 as observed in the RRL H53a with the VLA. The color 
scale shows the velocity field in the range from 5415 to 5615 km s^^. Contour levels are 
5425, 5445, 5605 km s~^ in steps of 20 km s~^. The crosses show the position of the two 
43 GHz continuum peaks. The synthesized beam is 0'.'7 FWHM. 



